Superposition of ferromagnetic and antiferromagnetic spin chains 
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Based on density functional theory band structure calculations, quantum Monte-Carlo simula- 
tions, and high-field magnetization measurements, we address the microscopic magnetic model of 
BaAg2Cu[V04]2 that was recently proposed as a spin-| anisotropic triangular lattice system. We 
show that the actual physics of this compound is determined by a peculiar superposition of fer- 
romagnetic and antiferromagnetic uniform spin chains with nearest-neighbor exchange couplings 
of ji^^ ~ — 19 K and Jp' ~ 9.5 K, respectively. The two chains featuring different types of the 
magnetic exchange perfectly mimic the specific heat of a triangular spin lattice, while leaving a 
clear imprint on the magnetization curve that is incompatible with the triangular-lattice model. 
Both ferromagnetic and antiferromagnetic spin chains run along the crystallographic a direction, 
and slightly differ in the mutual arrangement of the magnetic Cu04 plaquettes and non-magnetic 
VO4 tetrahedra. These subtle structural details are, therefore, crucial for the ferromagnetic or 
antiferromagnetic nature of the exchange couplings, and put forward the importance of comprehen- 
sive microscopic modeling for a proper understanding of quantum spin systems in transition-metal 
compounds. 

PACS numbers: 75.30.Et, 75.10.Pq, 71.20.Ps, 75.50.Gg 



I. INTRODUCTION 

Frustration and dimensionality are two crucial param- 
eters underlying the physics of magnetic systems. In in- 
sulators, these parameters rarely correlate with the ap- 
parent features of the atomic arrangement, because su- 
perexchange couplings are highly sensitive to details of 
the electronic structure and to positions of non-magnetic 
atoms linking the magnetic sites. While computational 
techniques based on electronic structure calculations de- 
veloped into a powerful tool for elucidating spin lattices 
of complex materials, simple phenomenological criteria 
are equally important for the preliminary assessment of 
the experimental data and the compound under consid- 
eration. 

The best-known phenomenological criterion of the 
magnetic frustration is the \6\/Tn ratio. It compares 
the Curie- Weiss temperature 9, which is often consid- 
ered as an effective energy scale of the magnetic cou- 
plings, to the magnetic ordering temperature T^fii High 
\6\/Tn ratios are believed to indicate strong frustration, 
although this rule will only hold for simple systems with 
few exchange couplings and well-established dimension- 
ality. Thus, the |6'|/rjv ~ 50 - 100 ratio is easily ob- 
tained even in non- frustrated quasi-one-dimensional (ID) 
systems, where strong quantum fluctuations due to the 
weak interchain couplings effectively prevent the system 
from long-range ordering down to low temperaturesi^^— 
Another possible scenario is that of magnets with strong 
dimer correlations, where the long-range-ordered state 
competes with the disordered singlet ground state, and 



the ordering temperature T^r may be strongly reduced 
without any frustration involvedi^^— The low \9\/Tiq ra- 
tio can be equally deceptive, because 9 is in fact a lin- 
ear combination of different exchange couplings that can 
be much smaller than the effective energy scale of the 
system. For example, the coexistence of ferromagnetic 
(FM) and antiferromagnetic (AFM) couplings renders 9 
and \9\/Tm low even in strongly frustrated magnets.— 

The phenomenological assessment of the frustration in 
a magnetic system has to be backed by additional cri- 
teria. Magnetic specific heat is an especially appeal- 
ing quantity, because it is expressed in absolute units 
and does not require an ambiguous reference to the ef- 
fective energy scale of the system. Further, the mag- 
netic specific heat distinguishes between the effects of 
dimensionality and frustration, with the latter leading 
to a much stronger reduction in the maximum of the 
magnetic specific heat {Cm)- For example, the spin-i 
square lattice (two-dimensional, non-frustrated) reveals 
the maximum of Cm/ R — 0.44, the spin-i uniform chain 
(one-dimensional, non-frustrated) shows a lower maxi- 
mum of Cm/ R — 0.35, but the specific heat maximum 
for the spin-i triangular lattice (two-dimensional, frus- 
trated) is even lower, Cm/R — 0.22iii The reduced mag- 
netic specific heat is a seemingly unambiguous measure 
of the frustration. It can be equally used to identify 
strongly frustrated spin systema^ '^^d^ or to refute pre- 
mature conclusions on the strong frustration^^ However, 
this phenomenological criterion is not universal, as we 
demonstrate below. 

In our study, the breakdown of the simple relation- 
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FIG. 1. (Color online) Top panel: perspective view of the 
BaAg2Cu[V04]2 structure showing alternating layers consist- 
ing of V04-bridged chains of Cul and Cu2, respectively. Dif- 
ferent colors (shadings) identify the inequivalent Cu04 pla- 
quettes and their slightly different orientation. Jc refers to 
the interlayer coupling. Bottom panel: a single layer in the 
ab plane (left) and the respective spin lattice with the intra- 
chain coupling Ja as well as interchain couplings Jati and Jab2 
(right). 



ship between the magnetic specific heat and the frustra- 
tion is related to a peculiar superexchange scenario in 
BaAg2Cu[V04]2. This compound has a fairly complex 
crystal structure with magnetic Cu^+ ions interspersed 
between the non-magnetic [V04]'^~ tetrahedra as well 
as Ba^+ and Ag+ cationsii^ The spatial arrangement of 
Cu^+ (Fig. [U resembles a weakly anisotropic triangular 
lattice with the intraplane Cu-Cu distances of 5.45 A 
(Ja), 5.63 A (Jabi), and 5.69 A {Jab2) and the interplane 
distance of 7.20 A ( Jc)- This lattice topology should in- 
duce magnetic frustration, as further corroborated by the 
magnetic specific heat that reaches the maximum value 
of Cm/R — 0.22 and strongly resembles theoretical pre- 
dictions for the spin-i triangular lattice 

In the following, we will show that the reduced Cm 
has a different origin, and arises from a peculiar super- 
position of FM and AFM spin chains. The system is, 
therefore, quasi-one-dimensional and only weakly frus- 
trated, in contrast to the straight-forward phenomeno- 
logical assessment. To support the one-dimensional sce- 
nario, we perform extensive band structure calculations 
combined with the fitting of magnetization and specific 
heat data. We also present original experimental results 
on the high-field magnetization that unequivocally rules 
out the triangular-lattice spin model for BaAg2Cu[V04]2. 



II. METHODS 

Our microscopic magnetic model of BaAg2Cu[V04]2 
is based on full-potential scalar-relativistic density func- 
tional theory (DFT) band structure calculations per- 
formed in the FPLO codoiS, implementing the basis set 
of local orbitals. We used the local density approxi- 
mation (LDA) with the Perdew-Wang parametrization 
for the exchange-correlation potentialJ^ The k meshes 
of 292 points and 64 points in the symmetry-irreducible 
part of the first Brillouin zone were chosen for the crys- 
tallographic unit cell and supercell, respectively. Corre- 
lation effects were treated on a model level or within the 
mean-field local spin-density approximation (LSDA)+t/ 
approach, as further described in Sec. IIIII 

Thermodynamic properties were calculated with the 
loopii and dirloop_sse (directed loop in stochastic se- 
ries expansion representation)^^ quantum Monte-Carlo 
(QMC) algorithms implemented in the ALPS simulation 
package Simulations were done for finite lattices with 
periodic boundary conditions. We used two independent 
chains containing L = 40 sites each. This chain length 
is sufficient to eliminate finite-size effects for thermody- 
namic properties within the temperature range under in- 
vestigation. 

Powder samples of BaAg2Cu[V04]2 were prepared ac- 
cording to the method described in Ref. [ij. Magnetic 
susceptibility was measured with MPMS SQUID magne- 
tometer in the temperature range 2 — 380 K in the ap- 
plied field of 0.1 T. Magnetization isotherm was collected 
at 1.5 K using the pulsed magnet installed in Dresden 
High Magnetic Field Laboratory. Details of the exper- 
imental procedure are described elsewherci^ The low- 
temperature heat capacity was measured above 0.5 K by 
a relaxation method using the '^He option of the Phys- 
ical Property Measurement System (PPMS, Quantum 
Design). 

III. MICROSCOPIC MAGNETIC MODEL 

LDA results for the band structure of BaAg2Cu[V04]2 
(Fig. [2]) closely follow expectations for a Cu^^-based insu- 
lating compoundi^"— Oxygen 2p states between —6 and 
—2 eV are surmounted by Ag 4c? and Cu 3d bands. The 
states above 2 eV originate from unfilled V id orbitals. 
While silver states are mostly found below —0.3 eV, Cu 
3d states additionally form narrow bands in the vicinity 
of the Fermi level. The calculated partial densities of 
states confirm the anticipated valences of Ag^+ (4^^"), 
Cu2+ (3d9), and (Sd"), and identify Cu2+ ions as the 
magnetic sites in the structure. The metallic LDA energy 
spectrum violates the insulating nature of the compound, 
as evidenced by the dark-yellow color of BaAg2Cu[V04]2. 
This discrepancy is well understood, given the impor- 
tance of correlation effects for the partially filled Cu id 
shell and the severe underestimation of such correlations 
in LDA. The missing correlations can be introduced on 
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FIG. 2. (Color online) LDA density of states 
BaAg2Cu[V04]2. The Fermi level is at zero energy. 



for 



the model level, or by a mean-field LSDA+C/ procedure. 

Following the first approach to the treatment of corre- 
lations, we consider in more detail the narrow bands in 
the vicinity of the Fermi level (Fig. [3]). The two bands 
can be assigned to two inequivalent Cu sites in the crystal 
structure. Both bands have the dx2^y2 orbital character, 
with X and y axes directed along shorter Cu-0 bonds. 
In BaAg2Cu[V04]2, the local environment of Cu^+ re- 
sembles a severely elongated octahedron, Cu04-)-2, with 
four short Cu-O bonds (1.96 — 1.97 A) lying in the plane 
and two long bonds (2.44 A) perpendicular to this plane. 
Therefore, the d^2_y2 orbital is the highest- lying crystal- 
field level in agreement with the LDA results. 

To fit the dx2-y2 bands with the tight-binding model, 
we construct Wannicr functions localized on Cu sites.— 
The fit perfectly reproduces the calculated band struc- 
ture (Fig. [3]), and yields Cu-Cu hopping parameters 
ti (Table |l|. Mapping the tight-binding model onto a 
one-orbital Hubbard model with the effective on-site 
Coulomb repulsion Uctf = 4.5 eVf^i^ we find the an- 
ticipated strongly correlated regime (ti <IC Ucs), and uti- 
lize second-order perturbation theory for analyzing the 
lowest-lying (magnetic) excitations. This way, AFM con- 
tributions to the exchange couplings are evaluated as 

The results of our model analysis are summarized in 
Table IP^ While AFM couplings in BaAg2Cu[V04]2 are 

(2) 

mostly weak, we find the sizable AFM coupling Ja along 
the a direction. Remarkably, this AFM coupling along 
a (denoted Ja) is observed for the Cu2 site and not for 
the Cul site, as emphasized by the superscripts (1) and 
(2) in the notation of Ji. This observation puts forward 
one important feature of the BaAg2Cu[V04]2 structure. 
The two Cu sites in BaAg2Cu[V04]2 are very similar and 
look nearly identical with respect to the geometry of in- 
dividual superexchange pathways (Table HI]). The Cul- 
Cul and Cu2-Cu2 distances in the ab plane are equal 
because of the constraints imposed by the lattice trans- 
lations. However, our microscopic analysis puts forward 




FIG. 3. (Color online) LDA bands (thin light lines) and the fit 
with the tight-binding model (thick dark lines). The k path 
is defined as follows: r(0,0,0), X(0.5,0, 0), M(0.5, 0.5, 0), 
F(0, 0.5,0), Z(0, 0,0.5), r(0.5, 0,0.5), i?(0.5, 0.5, 0.5), and 
yl(0, 0.5, 0.5), where the coordinates are given in units of the 
respective reciprocal lattice parameters. 



TABLE I. Cu-Cu distances (in A), hoppings ti (in meV), 
and exchange couplings Ji (in K) in BaAg2Cu[V04]2. The 
AFM contributions J:^™ 

are calculated as Ati/UeS with 
C/cff = 4.5 eV; the full exchange couplings Ji are obtained 
from LSDA+f/ calculations (Ud ~ 6 eV, Jd = 1 eV); and 
Jf^ = Ji- J^^^. The notation of J, is illustrated in Fig.[Tl 
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important differences between the deceptively similar su- 
perexchange pathways within the Cul and Cu2 sublat- 
tices (Table m. This difference gives a clue to understand 
the magnetism of BaAg2Cu[V04]2, and will be discussed 
in more detail below. 

The FM part of the superexchange originates from 
processes beyond the one-orbital model employed in our 
tight-binding analysis. In cuprates, FM interactions are 
generally ascribed to the Hund's coupling on the lig- 
and site2^ and can be evaluated by mapping total en- 
ergies for different coUincar spin configurations onto the 
classical Heisenberg model. The total energies arc ob- 
tained from spin-polarized band structure calculations 
with LSDA+[/ as the mean-field correction for correla- 
tion effects. Following previous studies of Cu^+-based 
compounds,— >2i we use the around-mean-field double- 
counting correction scheme, the on-site Coulomb repul- 
sion parameter Ud = 6 eV, and the Hund's exchange pa- 
rameter Jd = 1 eV. In the case of BaAg2Cu[V04]2, alter- 



4 




FIG. 4. (Color online) Wannier function based on the Cu 
d^2_y2 orbital. 



ations of Ud and the double-counting correction scheme 
have marginal influence on the results, and do not change 
the qualitative microscopic scenario. 

The total exchange couplings Ji based on the LSDA+J7 
calculations are listed in the last column of Table HI We 
find comparable FM contributions to the couplings ji^"* 

(2) 

and Ja along the a direction. Owing to the larger AFM 

(2) 

contribution to , this coupling remains weakly AFM, 
while ji^'' becomes FM. Other couplings show small FM 
contributions and hover around zero. The LSDA+JJ cal- 
culations confirm the leading couplings along a as well 
as the notable difference between Jq^' and Ja^\ Before 
comparing our magnetic model to the experimental data, 
we further comment on the microscopic origin of differ- 
ent exchange couplings in the Cul and Cu2 sublattices 
of BaAg2Cu[(V04]2- 

The sizable FM and AFM contributions are identified 
for the exchange couplings Ja^^ and Jp\ only. This find- 
ing is easily rationalized based on the magnetic dx2_y2 
orbital of the Cu^"'" ions. The crystal structure is best 
viewed in terms of the CUO4 plaquettes entailing the 
magnetic orbitals. This representation underscores the 
ID nature of the structure (Fig. [T]), and illustrates the 
quasi- ID magnetic behavior. However, unlike the well- 
known spin-chain Cu^+ compounds, such as Sr2Cu03 
(Ref.i) and CuPzN^i21 BaAg2Cu[V04]2 reveals a com- 
bination of two inequivalent spin chains with strikingly 
different exchange couplings. 

According to Tabled both Jq^' and Jq^' feature simi- 
lar FM contributions, yet very different AFM exchanges 
arising from different Cu-Cu hoppings in the effective 
one-orbital model. To elucidate the origin of these cou- 
plings, we consider the Wannier functions localized on 
Cu sites. Apart from the Cu d^2_y2 orbital forming the 
core of the Wannier function, we find sizable contribu- 
tions from oxygen 2p and vanadium 3c? orbitals (Fig.|4|). 
These contributions can also be observed in the LDA en- 
ergy spectrum (Fig. The Wannier functions of the 
neighboring Cu atoms overlap on the vanadium sites. 




FIG. 5. (Color online) Comparison of the Cul (left) and Cu2 
(right) chains in the BaAg2Cu[V04]2 structure. Note the dif- 
ferent orientations of the VO4 tetrahedra with respect to the 
Cu04 plaquettes, as quantified by the respective dihedral an- 
gles ip^^^ and (p^^\ 



TABLE II. Interatomic distances (in A) and angles (in deg) in 
the BaAg2Cu[V04]2 structure. The columns refer to the Cul 
and Cu2 layers, as shown in Fig.[T] The notation of individual 
atoms follows Fig. [S] (see text for details). 
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where each Wannier function features a different 3d or- 
bital of vanadium. This leads to the Hund's exchange on 
the vanadium site and explains the sizable FM contribu- 
tions to ji^' and in contrast to the very low FM 
contributions to other nearest-neighbor couplings hav- 
ing similar Cu-Cu distances (Table Note that a com- 
parable ~ — 15 K has been found in /3-CU2V2O7, 
where vanadium 3c? orbitals also contribute to the Cu- 
based Wannier functionsjSi 

We now consider different AFM contributions to ji^-* 
(2) 

and Ja ■ Geometrical parameters summarized in Ta- 
ble In] demonstrate a striking similarity between the re- 
spective superexchange pathways for Cul and Cu2. The 
only notable difference is the orientation of the VO4 
tetrahedra with respect to the chains. Naively, the po- 
sition of the tetrahedra is described by the 08-02-01 
and 07-04-03 angles (Fig. [5]). However, these do not 
account for the different tilting of the CUIO4 and Cu204 



plaquettes with respect to the a axis (Fig. [T|). Therefore, 
we use dihedral angles (p referring to the 01 '-08-02 and 
01-02-01'-02' planes for Cul {ip^^'>), and to the 03'- 
07-04 and 03-04-03'-04' planes for Cu2 
cording to Table HIl the difference between (^^^^ and (^^^^ 
is as large as 21.5 deg, thus to be considered as the main 
feature to account for the different Cu-Cu hoppings ii^-* 
and ii^''. 

To explore the role of the dihedral angles ip, we con- 
struct fictitious model structures with the VO4 tetrahe- 
dra rotated about the 0-0 edges (01-02 and 03-04 
for VI and V2, respectively). This way, we are able to 
tune 93*^^^ toward (/s*-^^ = 102.2 deg and enhance ii^-* to 
21 meV (compare to —11 meV at the experimental (/s*-^^ = 
123.7 deg), or change ip^-^^ toward (p^^^ = 123.7 deg, thus 

reducing ta to 3 meV (compare to —43 meV at the 
experimental v?^^^ = 102.2 deg). Overall, a change of ori- 
entation by approximately 22 deg is accompanied by a 
A{ta) of 32 meV and 46 meV, respectively. Therefore, 
the orientation of the non-magnetic VO4 tetrahedra is 
of crucial importance for the Cu-Cu hoppings and AFM 
superexchange. Note, however, that this geometrical pa- 
rameter is not unique, and the specific arrangement of 
the CUO4 plaquettes with respect to the chain direction 
(Fig. [T] and Fig. [S]) is also responsible for the large AFM 

contribution to Ja ■ 



IV. EXPERIMENTAL DATA 

The DFT results summarized in Table U identify the 
spin lattice of BaAg2Cu[V04]2 as a system of weakly 
interacting inequivalent spin chains with the intrachain 

couphngs Ja^^ and respectively. While J^^^ is clearly 
(2) 

FM, Ja is weakly AFM and probably close to zero. 
This qualitative scenario is verified by the magnetization 
isotherm measured at 1.5 K. Previous measurements^^ in 
fields up to 5 T showed that half of the Cu spins seem to 
saturate around 1.5 T. Here we extend our study into the 
behavior of the magnetization in higher fields (Fig. 
Based on these high-field measurements, we show that 
the magnetization of BaAg2Cu[V04]2 is further increased 
between 1.5 T and 16 T, where the full saturation with 
M ~ 1.08 /is/f.u. is reached. This peculiar behavior ap- 
parently contradicts the conjecture on the triangular spin 
lattice that would lead to a smooth increase in the mag- 
netization between zero field and the saturation fieldi^ 
The experimental magnetization curve is readily elu- 
cidated by our microscopic model. While half of the 
spins comprising the FM spin chains (Cul) align with 
the field already at 1.0 — 1.5 T once thermal fluctuations 
are suppressed, the remaining spins (Cu2) are coupled 
antiferromagnetically and require larger fields to over- 
come the AFM interactions. This behavior strongly re- 
minds of a two-sublattice ferrimagnet, where half of the 
maximum magnetization is recovered in low fields, while 
larger fields arc required to flip one of the sublattices. 



5 



^ 0.4 





' 1 ' 
Fit (FM+ AFM chains) 


1 ' 


- 


Experiment (r = 1 .5 K) 








FM cliain, Jj" = -19K 






AFM ciiain, jj'* = 9.5 K 



10 20 30 

Field (T) 

FIG. 6. (Color online) Magnetization isotherm of 
BaAg2Cu[V04]2 measured at 1.5 K (filled circles) and the fit 
with a combination of FM and AFM spin chains (solid line). 
The contributions of the FM (Cul) and AFM (Cu2) chains 
are shown by the dashed and dotted lines, respectively. 
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FIG. 7. (Color online) Magnetic susceptibility of 
BaAg2Cu[V04]2 measured in the applied field of 0.1 T (filled 
circles) and the fit with a combination of FM and AFM spin 
chains (solid line). The contributions of the FM (Cul) and 
AFM (Cu2) chains are shown by the dashed and dotted lines, 
respectively. 

Note, however, that BaAg2Cu[V04]2 is not in a mag- 
netically ordered state at 1.5 K, hence no magnetization 
hysteresis is observed. The long-range magnetic order in 
BaAg2Cu[V04]2 is estabhshed below Tc ~ 0.7 K and is 
further discussed in Sec. |Vl 

The above qualitative picture can be quantified 
by fitting the experimental magnetization data.— In 
BaAg2Cu[V04]2, field dependence of the magnetization 
(Fig.ini) and temperature dependence of the susceptibility 

(Fig. [T]) are complimentary. The magnetization isotherm 

(2) 

is sensitive to the AFM exchange Js, that determines the 
saturation field, but the alignment of the FM component 
mostly depends on thermal fiuctuations so that Ja^-* can- 
not be determined precisely. In contrast, the FM chains 
coupled by J^^^ produce the dominant contribution to the 
susceptibility;^ which gives an accurate estimate for Jq^-* , 
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FIG. 8. (Color online) Magnetic part of the specific heat 
[Cm{T)] divided by the gas constant (R) for BaAg2Cu[V04]2 
measured in zero field (top) and in applied fields of 3 T (mid- 
dle) and 7 T (bottom). The simulated curves for the com- 
bination of FM and AFM spin chains are shown by solid 
lines, whereas the dashed and dotted lines denote the con- 
tributions of the FM (Cul) and AFM (Cu2) spin chains, re- 
spectively. Experimental data (circles) are taken from Ref.[l3. 
The model parameters ji^' = — 19 K and jI^' = 9.5 K are 
extracted from the fits to the magnetization data (Figs. |6] 
andO. Therefore, we compare our model to the experiment 
with no adjustable parameters. 



tibility of BaAg2Cu[V04]2, respectively. The FM chains 
lead to the sharp increase in the susceptibility at low tem- 
peratures, while the contribution of the AFM chains is 
barely visible on the same scale. The contribution of the 
FM chains to the magnetization isotherm is saturated at 
low fields and corresponds to one half of the maximum 
magnetization, because half of the Cu atoms belong to 
the FM chains. The magnetization of the AFM chains is 
linear at low fields, bends upward above 7 T, and finally 
saturates around 16 T where the full alignment of spins 
is achieved. 

Wc will now test our quasi-lD model against the exper- 
imental specific heat data showing the strongly reduced 
maximum that might be characteristic of a spin-i trian- 
gular lattice. We use the fitted parameters based on the 
magnetization data and, therefore, compare our model to 
the experiment with no adjustable parametersi^ Fig. [5] 
presents the magnetic specific heat data measured in zero 
field and in two representative applied fields along with 
the simulated curves. The remarkable agreement be- 
tween the experiment and the model prediction confirms 
our microscopic scenario, and suggests that the strongly 
reduced specific heat maximum, especially in zero field, 
is not an unambiguous footprint of the magnetic frustra- 
tion. 

In zero field, the specific heat maximum closely follows 
the contribution of the AFM spin chains, while the FM 
chains with the stronger coupling Ja^^ ~ — 19 K provide 
a temperature-independent "background" below 15 K. 
The applied field of 3 T increases the maximum up to 
Cjn/R ^ 0.32. The stronger field of 7 T additionally 
shifts the maximum to higher temperatures. Both ef- 
fects are perfectly reproduced by our microscopic model. 
Magnetic fields transform the temperature-independent 
zero-field specific heat of the FM chains into a small max- 
imum at 3.5 — 4.0 K. This maximum of the FM contribu- 
tion weakly depends on the field, because the FM (Cul) 
subsystem is saturated above 2 T (Fig. |6]). By contrast, 
the contribution of the AFM chains shows a pronounced 
field dependence that underlies the evolution of the ex- 
perimental magnetic specific heat in the applied field. 



while leaving certain ambiguity for Ja ■ The two sets of 
data are successfully fitted with the same model param- 
eters: ~ -19 K, ^ ~ 9.5 K, g ~ 2.16 (Figs. [6] 
and [7]). We also included a temperature-independent 
contribution to the susceptibility xo — 9 x 10^^ emu/mol 
that accounts for the van Vleck paramagnetism and core 
diamagnctism. Our fitted g-value is in excellent agree- 
ment with the experimental powder-averaged g = 2.18f^ 
While J^^^ closely follows the DFT prediction (Table IJ), 

(2) 

the computational estimate of Ja is less accurate, al- 
though still acceptable considering the low energy scale 
of the exchange couplings in BaAg2Cu[V04]2'^ 

Figures [H] and [7] illustrate the contributions of the FM 
and AFM components to the magnetization and susccp- 



V. DISCUSSION AND SUMMARY 

The combination of DFT calculations and QMC fits to 
the experimental data gives compelling evidence for the 
quasi-lD magnetic behavior of BaAg2Cu[V04]2. The su- 
perposition of FM and AFM spin chains with different 
magnitudes of the exchange couplings results in pecu- 
liar and perplexing thermodynamic properties. While the 
zero-field specific heat resembles the typical response of 
the spin-i triangular lattice, the magnetization isotherm 
is reminiscent of a system with two different magnetic 
sublattices, and underpins the proposed magnetic model. 

Based on our microscopic analysis, we establish the 
spin lattice of BaAg2Cu[V04]2 as a peculiar derivative 
of conventional Heiscnberg spin chains with nearest- 
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FIG. 9. (Color online) Field-dependent magnetic part of the 
specific fieat, Cm{T), divided by R for BaAg2Cu[V04]2 at 
low temperatures, showing the behavior typical for a ferro- 
magnetic or ferrimagnetic transition. 



neighbor exchange coupling J. This model was widely 
studied for both FM and AFM J^'^^ but the combina- 
tion of FM and AFM spin chains was neither considered 
theoretically nor encountered experimentally. 

The superposition of inequivalent spin chains is a chal- 
lenge for "non-local" experimental techniques, such as 
thermodynamic measurements or inelastic neutron scat- 
tering, that probe the system as a whole. These methods 
inevitably blend the signals of different sublattices, and 
generally lead to a complex response that can be fully 
elucidated based on the microscopic approach only. A 
more direct experimental information could be extracted 
from "local" methods, which probe different magnetic 
sublattices independently. For example, an elegant way 
to study the physics of BaAg2Cu[V04]2 further could be 
nuclear magnetic resonance (NMR) on ^^V nuclei. The 
inequivalent vanadium sites VI and V2 are coupled to 
Cul and Cu2, respectively. Owing to the very similar 
local environment, the signals from these two vanadium 
sites should perfectly overlap at high temperatures. At 
low temperatures though, the lines will split because of 
the different Knight shifts resulting from the disparate 
local magnetization in the vicinity of the FM and AFM 
spin chains. Therefore, the NMR experiment should be a 
valuable additional experimental test of our microscopic 
model. 

Another interesting problem is the long-range-ordered 
(LRO) ground state of BaAg2Cu[V04]2. While isolated 
spin chains do not show the LRO even at zero tempera- 
ture, interchain couplings induce the LRO state at a finite 
temperature)^"— irrespective of the weak frustration 
that could be induced by the triangular arrangement.— 
In BaAg2Cu[V04]2, specific- heat measurements reveal 
the sharp anomaly at Tq — 0.7 K in zero field. This 
anomaly is drastically suppressed even in weak magnetic 



fields (Fig. |9l see also Ref. [13), as typical for a ferromag- 
netic transition, or - more generally ~ for an LRO state 
with non-zero net magnetization. Such ground state can 
be indeed derived from our microscopic model and ex- 
plained in terms of a two-sublayer system with the inter- 
layer exchange coupling Jc- As outlined above, these sub- 
layers are stacked along the c-axis in an alternate fashion. 
Each plane consists either of FM (Cul) or AFM (Cu2) 
spin chains, respectively (Fig. [T]). 

The Cul spins within the FM spin chains prefer the 
parallel alignment so that a FM sublattice is formed. The 
Cu2 spins are expected to be ordered antiferromagneti- 
cally along a and form an AFM sublattice. The nature 
of the interchain couplings is more difficult to establish 
because of their lower energy scale that might allow for 
additional, non-isotropic contributions, such as dipolar 
interactions. However, even the isotropic (Heisenberg) 
model based on DFT enables us to make a plausible con- 
jecture about the ground state. The FM couplings Jabi 
and Jab2 in the ab plane are compatible with both FM 
and AFM exchange along a. These couplings should re- 
inforce the formation of the FM sublattice for Cul and 
the AFM sublattice for Cu2. The AFM coupling Jc along 
c introduces a weak frustration, but its effect should be 
small. In summary, BaAg2Cu[V04]2 entails two inequiv- 
alent sublattices and presents a peculiar example of a 
spin-i system with non-zero net magnetization. 

From phenomenological point of view, a similar ground 
state with the non-zero net magnetization has been re- 
cently observed in the spin-i fcrrimagnet Cu20Se034S 
However, unlike conventional ferrimagnets and unlike 
Cu20Se03, BaAg2Cu[V04]2 does not feature well- 
defined sublattices with opposite directions of the spin, 
and rather shows a sequence of FM and AFM layers. Un- 
fortunately, the frustrated nature of the interlayer cou- 
pling Jc prevents us from using QMC for simulating the 
ground-state properties and the transition temperature 
Tc- Therefore, we are presently unable to verify the pro- 
posed magnetic structure. Further experimental studies, 
such as neutron diffraction, would be required to tackle 
this problem. 

The microscopic magnetic model of BaAg2Cu[V04]2 is 
furthermore instructive from a structural viewpoint. The 
Cul and Cu2 sites look deceptively similar, so that one 
would not expect any substantial difference between the 
magnetic couplings within the two sublattices. However, 
the couplings are very different - not only in the magni- 
tude but also in the nature - because of the subtle infiu- 
ence of the VO4 tetrahedra connecting the neighboring 
CUO4 plaquettcs. The effect of the non-magnetic group is 
sizable and two-fold. Vanadium 3c? orbitals contribute to 
the Wannier functions, and induce a FM superexchange, 
which is weakly dependent on the specific arrangement 
of the VO4 tetrahedra. This FM contribution represents 
a constant term that is superimposed on a variable AFM 
superexchange. The latter is controlled by the Cu~Cu 
hoppings showing dramatic dependence on the mutual 
orientation of the VO4 tetrahedra and CUO4 plaquettes. 
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Depending on the specific geometry, the AFM contribu- 
tions may or may not surpass the FM superexchange, 
and quahtatively different exchange couplings emerge. 

The subtle dependence of AFM superexchange on the 
orientation of non-magnetic tetrahedra is reminiscent of 
the long-range couplings in BiCu2P06, where slight ro- 
tations of the bridging PO4 groups modify the interac- 
tions by 50 — 70 K^i More generally, the unusual micro- 
scopic scenario of BaAg2Cu[V04]2 confirms the crucial 
importance of non-magnetic bridging groups for the su- 
perexchange in magnetic insulators. Other remarkable 
examples include the effect of Ge04 tetrahedra on the 
Cu-based spin chains in CuGeOa (Ref. IH) , as well as the 
unusual ferromagnetism of CdVOs related to the low- 
lying 5s orbitals of Cd atoms.— The effect of the non- 
magnetic groups opens broad prospects for tweaking su- 
perexchange couplings by minor alterations of the crystal 
structure. For example, BaAg2Cu[V04]2 is likely to sus- 
tain cation substitutions in the Ba and Ag positions, thus 
leading to further interesting combinations of FM and / or 
AFM spin chains in a single chemical compound. 

In summary, we have derived a microscopic magnetic 
model of BaAg2Cu[V04]2, and presented a consistent in- 
terpretation of the available experimental data for this 
compound. The crucial and highly unexpected feature of 
BaAg2Cu[V04]2 is the dramatic difference between the 
couplings within the Cul and Cu2 sublattices. While the 
Cul sublattice is ferromagnetic, the Cu2 sublattice is an- 



tiferromagnetic. This unusual - and so far unreported - 
combination of weakly coupled FM and AFM spin chains 
within a single chemical compound leads to peculiar ther- 
modynamic properties, with the specific heat resembling 
that of a strongly frustrated two-dimensional spin sys- 
tem. The spin lattice of BaAg2Cu[V04]2 is, however, 
only weakly frustrated and quasi- ID, as confirmed by the 
high-field magnetization measurements suggesting the 
ground state with non-zero net magnetization. The dif- 
ferent couplings within similar structural units are solely 
determined by the orientation of the non-magnetic VO4 
tetrahedra with respect to the CUO4 plaquettes. These 
results present an instructive example on the importance 
of bridging groups for superexchange pathways, and open 
interesting opportunities for tuning low-dimensional spin 
systems within a given structure type. 
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